Differential Renormalization-Group Approach 
to the Layered sine-Gordon Model 

I. NandorJ] 

Institute of Nuclear Research of the Hungarian Academy of Sciences, 
H-4001 Debrecen, P.O.Box 51 Hungary; 
Max Planck Institute for Nuclear Physics, 

in ! 

O- P.O.Box 103980, Heidelberg 69029, Germany 

O 1 

<n : 

OJO' K. Sailer 

3 : 

^ ■ Department of Theoretical Physics, University of Debrecen, 

H-4010 Debrecen, P.O.Box 5 Hungary 
(Dated: February 7, 2008) 

cn ■ 

cn ! Abstract 
O ■ 

New qualitative picture of vortex length-scale dependence has been found in recent electrical 
transport measurements performed on strongly anisotropic B^S^CaC^Os single crystals in zero 
magnetic field. This indicates the need for a better description of the 3D/2D crossover in vortex 
Qh! dimensionality. The vortex-dominated properties of high transition temperature superconductors 

with extremely high anisotropy (layered systems) are reasonably well described in the framework of 



X 



the layered XY model which can be mapped onto the layered sine-Gordon model. For the latter we 
derive an exact renormalization group (RG) equation using Wegner's and Houghton's approach in 
the local potential approximation. The agreement of the UV scaling laws find by us by linearizing 
the RG equations with those obtained previously in the literature in the dilute gas approximation 
makes the improvement appearant which can be achieved by solving our RG equations numerically. 
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I. INTRODUCTION 



Vortex dominated properties of high-temperature superconductors (HTSC) are usually 
investigated by means of electrical transport measurements which have raised many ques- 
tions about the critical behaviour and the length-scale dependence of vortices in layered 
systems. The HTSC materials have a layered structure and usually consist of copper-oxide 
superconducting planes separated by insulators. In case of extremely high anisotropy like 
in Bi 2 Sr 2 CaCu208 (BSCCO) single crystals, the coupling between the Cu02 planes is weak. 
Due to this weak coupli ng, one of the commonly accepted models is the layered XY model 



(Pierson and Vail 



199 2]) where the weak interlayer coupling is given by a Lawrence- Doniach 



type term ([Pierson et. al 



1992J). The topological excitations in superconducting layers are 



vortex-antivortex pairs. Two such pairs belonging to neighbouring layers can form vortex 
loops due to the Josephson coupling. 



(Pierson and Vails 


1994 


model ( 


Pierson 


1995 





1992) or it can be transformed into the layered Coulomb gas (LCG) 



1995). The critical behaviour of vortices and vortex loops in the LSG and 



LCG models have b een investigated by renormalization group (RG) methods using the dilute 



gas approximation (Pi erson 



1995 



Pierson et. al 



19921 ) which assumes low vortex fugacity. 



However, the vortex fugacity becomes large above the critical temperature (but below the 
Ginzburg-Landau transition temperature) when individual vortices can be spontaneously 
created although the sample is still in the s uperconducting p hase. 



Recent current transport measurements ()Vad et al 



20041 1 performed on BSCCO super- 



conducting single cry stals show a new qu alitative picture of vortex behaviour above the 



critical temperature (|Nandori et. al 



20041 ) . Using a multi-contact 'DC flux transformer' 



configuration, the temperature dependence of primary and secondary voltages are measured 
in zero magnetic field. (The side of the crystal where the current is injected is called the 
primary one, its opposite side is the seco ndary one.) Ac cording to the theoretical predictions 
based on the dilute gas approximation (jPiers on 1995), above the critical temperature one 
should expect the decrease of the secondary voltage to zero as a function of t he temperature. 
Instead, the experimental data reveal a new feature ((Nandori et. al\ 120041 ) . a second peak 
appears in the temperature dependence of the secondary voltage. A possible reason for the 
disagreement between the theory and experiment is the usage of the dilute gas approxi- 
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mation. This indicates the need for a better description of length-scale dependence of the 
vortex dimensionality at the 3D/2D crossover and hence the modification of the RG analysis 



(Nandori et. al 



2004). 



The purpose of this paper is to propose a theoretical study of t he LSG model using 



the differential RG approach ( Wilson and Kogut 



1974 



Wilson 



19831 1 in order to improve 



the treatment of the length-scale dependence of vortices. In section 2, the definition of 
the layered models and their interrelati ons are summarized. After gi ving a short sketch of 
Wegner's and Houghton's RG method ((Wegner and Houghton! 119731 1 for models with two 
interacting scal ar fields in sec tion 3, we apply it to the LSG model using the local potential 



approximation fPo lonvi 



2001) in section 4. The UV scaling laws are found and discussed in 



section 5 on the base of the linearized RG equations. The comparison of the UV sca ling laws 



to those obtained previously in Refs. ((Pierson and Vail 



199 



Pierson et. al 



1992) using a 



momentum-space RG method with smooth cut-off in the dilute gas approximation is also 
given, and the possibility of an improved treatment of the RG flow by our method is pointed 
out. Finally the results are summarized in section 6. 

II. LAYERED MODELS 



The phase structure of HTSC materials is usually investigated in the framework of the 
Ginzburg-Landau (GL) theory of superconductivity at the phenomenological level. The GL 
theory is based on a variational method which has been applied to the gradient expansion of 
the free-energy in powers of the complex order parameter ip and its derivatives where \ip\ 2 
represents the local density of superconducting electron pairs. In case of strong an isotropy 



and in the absence of external fields the free energy becomes ((Nandori et. al 



2004J) 



F = Y J j Q^||(V0 n ) 2 + J±[l ~ cos(0 n - ^OIJ 



(1) 



where ip n (r) = ip exp[i<p n (r)} is assumed with ip real and identical in every layer n. The 
coupling constants of the model are Jn = lftil) 2 /m a b and J± = h 2 tfj 2 /m c s 2 with the effective 
masses m a b,m c , and the interlayer distance s. The last term is related to the Lawrence- 
Doniach term, in which the coupling between the layers is the Josephson coupling. Expand- 
ing the cosine term in Taylor series around zero, one arrives at the layered two-dimensional 



3 



XY model (jPierson and Vallslll992h : 

Hlxy = \ £ / d 2 r J|,(V0„) 2 + \j2j d2rJ ^ 



>n-\) 



(2) 



The l ayered XY model c an be mapped onto the LSG or LCG models ()Pierson and Vails 



1994. 



Pierson et 



model reads as 



J LSG 



al 



1992) . For a system consisting of two layers the Lagrangian of the LSG 



Pierson and Vails 



19921 ) 



-[(d^) 2 + («90 2 ) 2 ] + - 2 ) 2 + u [cos(/fyi) + cos(/30 2 )] 



(3) 



where 0i, 2 are one-component scalar fields, u is the fugacity parameter for the vortices of 
the model, (3 2 = Jii and J = J±/(A 2 J\\) are dimensionless coupling constants and A stands 
for the UV momentum cutoff. The second term on the r. h.s. d escribes the interaction of 

19921 ) RG flow equations, 



the layers. In Refs. ((Pierson and Vail 
1 



1992 



Pierson et. al. 



k^-p 2 = — ( u 1 ^ - 4.5(3% Z J In J), 
dk 47r 



k±J=-2J, (4) 



where k is the momentum scale, have been obtained for the coupling constants of the LSG 
model (HH) with the help of a momentum space RG analysis using a smooth (anisotropic) 
momentum cutoff and the dilute gas approximation. The latter assumes that the fugacity 
of the vortex gas, i.e. the coupling constant m in Q is small. 

III. DIFFERENTIAL RG APPROACH 



The renormalization of the LSG model is presented in this section b y means of the differ- 



ential RG approach perfo rmed in momentum space with s 



blocking transformations (W ilson and Ko gut 



mi 



Wilson 



harp 



cutoff ([Polonvi 



200l|). The 



19831 ) are realized by successive 



elimination of the field fluctuations according to their decreasi ng momentum in infinitesima l 
steps. We use Wegner's and Houghton's method (WH-RG) ([Wegner and Houghto n 1973) 
when the high-frequency modes are integrated out above the moving sharp momentum cut- 
off k and the physical effects of the eliminated modes are built in the scale-dependence of 
the coupling constants. 

Assuming that the generating functional does not change under the infinitesimal blocking 
transformation when the momenta of the eliminated modes are taken from a thin momentum 
shell with radius k and thickness Ak, the blocking relation 



V\ 



-S k {<j>+4>'] 



(5) 



is obtained for the blocked action Sk[4>] (|Polonvill200lt ). The field variables and 0' contain 
Fourier modes with momenta p < k — Ak and k — Ak < p < k, respectively. Expanding the 



action [4> + 4>'} in Taylo r series around its sadd 



and the WH-RG method ((Wegner and Houghton 



e poin t, the path integral can be evaluated 



19731 ) provides the functional RG equation 



S k - Ak [4>} = S k [4>] + -trln(S'M) + 0(h 2 ) 



(6) 



where S'l denotes the second functional derivative of the blocked action with respect to the 
field and the trace is over the momentum shell [k — Ak, k]. In the limit Ak — > Eq. ((^J) 
reduces to an exact integro-differential equation for the functional Sk <j)} because the small 
parameter Ak/k suppresses the higher- loop contributions ( PolonvilboOl ). The generalization 



of Eq.floT) for two, interacti ng scalar fields is straightforward following the method proposed 



in (|Liao and Polonvi 



1995T ) and results in 



where S%[ 
to d>j and 



denote the second functional derivatives of the blocked action with respect 



3 ' 



If the saddle point is non-trivial in (jSJ), a spinodal instability occurs at some critical 
mom entum scale k c , and the WH-RG equation (JHJ) loses i ts validity ([A lexan dre and Polonvi 
1999 ). Then, the blocking produces a tree-level effect ( Alexandre and Polonvi T1999L In 
case of two interacting scalar fields the tree-level blocking relation reads as 



S k -Ak[<Pi,H= min (5 fc [0i + 0f,0 2 + ^]) 

±cl J,cl 



(8) 



where the minimum is sought over field configurations including only Fourier-modes with 
momenta from the infinitesimal momentum shell k — Ak < p < k. 



IV. RG APPROACH TO THE GENERALIZED LSG MODEL 

In this section we specify the WH-RG equation (J7J) and the tree-level blocking relation (jSJ) 
for the LSG model. First, the WH-RG equation and the tree-level blocking relation, have 
to be projected into a particular functional subspace, in order to reduce the search for a 
functional (the blocked action) to that of functions. Here we assume that the blocked action 
contains only local interactions. Expanding it in powers of the gradients of the fields <pi and 
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2 and keeping only t he leading order terms one arrives at the local potential approximation 



(LPA) (jPolonvi 



2001). In LPA the ansatz for the blocked action for two scalar fields can be 



written as 



d d x 



with the blocked potential V k 

kdkV k (<l>i, <h) = -k d a d \n 



2(dM 2 + -(dM 2 +V k { 



(9) 



2 ) satisfying the WH-RG equation 

fc 4 



(10) 



±0,(2 



where fid denotes the solid angle in d dimensions and Vl J 



With a d = pL d {-Z1T) , vvxxcxc *u d UCUUra me duuu cuxxgxc m u, muiaioiuup o,^ r fe 

d ( j >i d ( j > ,Vk{4>iT4>2)- I n order to remove the trivial scale-dependence of the coupling constants, 
the WH-RG equation (jlOj) has to be rewritten in terms of dimensionless quantities. For 
dimensions d = 2, the dimensionless form of Eq. (fTU|) reads as 



[2 + kd k )V k 



a 2 ln([l + 0[l + V; 22 ]-^ 21 



(11) 



where V k {4>i, <j>i) = k~ 2 Vk(4>±, ^2) h as been introduced. 

The argument of the logarithm in Eqs. (|10j) and ()lip must be positive. If the argument 
vanishes or if it changes sign at a critical value k c , the WH equation ([TU|) loses its validity 
for fc < k c . This is a consequence of the spinodal instability, i.e. that of the vanishing 
of the restoring force for field fluctuations of momenta from the infinitesimal momentum 
shell. Such fluctuations can grow to those of finite amplitude and result in a non-vanishing 
sadlle-point configuration. Below the critical scale k < k c the tree-level blocking relation (JHJ) 
can be used to find the saddle point. Looking for it among the configurations 4>f and <pf 
having the form of plane waves one finds for the blocked local potential 



Vi 



k-AkWl, 



mm 

pk 



2 Pk kZ + 7; I dp\40i + 2p k cos(ttp), <fi 2 + 2p k cos(np)) 



-1 



;i2) 



where p k is the amplitude of the plane waves ((Alexandre and Polonvilll999l ). Notice that the 
layers in the LSG model are physically equivalent therefore the potential is invariant under 
the transformation <pi <-> (p2 which implies identical plane-wave forms of (f>f and cf>2 with the 
same amplitude p k and wave vector in the saddle point. For dimensions d = 2, the tree-level 
relation (|T2*j) for dimensionless quantities reads as 

-1 



fc-AHrl> 



mm 

pk 



2pl + 



dpV r fc (0i + 2p k cos(7rp), 4> 2 + 2p fc cos(7rp)) 



-1 



(13) 
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As to the next, we specify the WH-RG equation (|TT|) and the tree-level blocking relation 
(fT5Jl for the LSG model. On the one hand, the ansatz for the blocked potential should 
be reach enough in order to ensure that the RG flow does not leave the chosen subspace 
of blocked potentials. On the other hand, the ansatz for the blocked potential should 
preserve all symmetries of the original model at the UV cutoff scale k = A. In particular, 
the blocked potential for the LSG model should be invariant under the exchange of the 
field variables, 0i <-> 02 because the layers are physically equivalent and also reveal the 
symmetries fa ~^ — 0i an d fa — > 4>i + jp- (c.f. Eq. (J3J)). Therefore, we choose the ansatz 

V k (fa, fa) = \Ufa - fa) 2 + U k (<f>i, fa), (14) 
where Uk(4>i, fa) is an arbitrary periodic function of the fields 

£4(0i, fa) = ^ [u nm {k) cos(n/30i) cos(m/50 2 ) + v nm (k) sm(n(3fa) sm(m(3fa)] (15) 

n,m 

with symmetric Fourier amplitudes u nm = u mn , v nm = v mn . The potential (fTHjl defines the 
generalized LSG model which contains higher harmonics of the field variables. Notice, that 
the scale-dependence is entirely encoded in the coupling constants Jk,u nm (k) and v nm {k). 
Inserting this ansatz into the WH-RG equation (fTTj) . 

(2 + kd k ) Q J k (fa - fa) 2 + U k (fa, fa)j = (16) 
-a 2 In (l + 2 J k + ill' + Uf + J h (U? + Uf + UD + U?V? - (C/ fc 12 ) 2 ) , 

and separating the periodic and non-periodic parts, one arrives at the RG flow equations 

(2 + kd k ) J k = 0, (17) 
(2 + kd k ) U k (fa, fa) = -a 2 In (l + 2 J k + U 1 , 1 + Uf + J k {Uf + Uf + U 12 ) 

+ ul l ul 2 - {ul 2 ) 2 ) 

since the argument of the logarithm is a periodic function of the fields. Notice that the 
dimensionless coupling constant J k has only a trivial tree-level evolution and hence the 
dimensionful J = k 2 J k remains constant during the RG procedure. It is also important to 
note that the second equation of (j!7j) keeps the periodicity with the same length of period, 
therefore j3 has no scale-dependence in the LPA. 
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In order to specify the tree-level blocking for the generalized LSG model, one has to 
insert the ansatz ()14j) and (JT5|l into the tree- level relation (jl3j) . For the sake of simplicity 
we turn back to dimensionful parameters and the relation (|12|) . Performing the integral and 
separating the periodic and the non-periodic parts, Eq. (J 12)) reduces to two relations. One 
of them, 



Jk- 



AA: 



Jk 



(18) 



says that the dimensionful coupling constant J remains unchanged under the tree-level 
blocking, so that the corresponding dimensionless coupling keeps its trivial scaling J k = Jk~ 2 
in LPA even below the critical scale k c . The relation for the periodic part is 



£4-Afc(0i, <h) = mm {2p 2 k {k 2 + J) + (19) 
pk 

22 [ u nm{k) cos(n/30i) cos(mf3(j)2)Jo(2nf3p k )J (2mf3p k ) 

n,m 

+ v nm (k) sin(n/30i) sin(m/?0 2 )^o(2n ( 5pfc) J (2m/?p fc )] } , 

where Jq(x) stands for the Bessel function. 

Eqs. (|17J) and (|T9"j) can be used for the determination of the RG flow of the generalized 
LSG model in the LPA, but they are solvable only numerically. Instead, we give below an 
analytic solution of the linearized RG equations which is valid only in the UV scaling regime. 
Nevertheless, it has the value of being closely related to the re sult obtained in the dilute 



gas approximation ((Pierson and V ails 



1992 



Pierson et. al. 



1992) and with its help one can 



recognize the improvement to be expected form the numerical determination of the WH-RG 
flow in the LPA. 



V. LINEARIZED RG 

In order to discuss the linearized RG flow of the generalized LSG model, the dimensionless 
WH-RG equation fTTf obtained in LPA is linearized around the Gaussian fixed point 

(2 + kd k ) V k (<f) U 2 ) = - a 2 (V k u + Vf ) (20) 



S 



by assuming 1 >> \V k lJ \ where V£ 3 = d^d^ Vfc(0i, 2 ). Using the ansatz (jT4*)l and (|T5|) for 
the generalized LSG model, the derivatives of the potential are 

V k u = J k - n 2 /5 2 [« nm (fc) cos(n/50x) cos(m/30 2 ) + v nm (k) sin(n/30i) sin(m/30 2 )] , (21) 

^jfc 22 = -4 - ^ ™ 2 /3 2 [«nm(^) cos(n/50i) cos(m/30 2 ) + sin(n/30i) sin(m/50 2 )] . 

Inserting Eqs. (J2T|) into the linearized WH-RG Eq. (f^Oj) the following RG equations are 
obtained for the coupling constants of the generalized LSG model 

(2 + kd k )J k = 0, (22) 
(2 + kd k )u nm {k) = a 2 (3 2 {n 2 + m 2 )u nm (k), 
(2 + kd k )v nm (k) = a 2 (3 2 (n 2 + m 2 )v nm (k), 

where a 2 = l/Air. The linearized RG equations for the various coupling constants decouple 
and their solutions can be obtained analytically 

J k = J k =A (r^j ' (23) 

i\ -2+a 2 f3 2 (n 2 +m 2 ) 

u nm [k) = u nm {k = A) 



Vnmik-) — V nm {k — A) 



-2+a 2 /3 2 (n 2 +m 2 ) 



A 
k_ 

where J k =A,u nrn (k = A) and v nm (k = A) are the initial values for the coupling constants at 
the UV cutoff A. The linearized RG flow predicts two phases of the model. If the parameter 
(5 2 is larger than a critical value j3 2 > /3 2 = 8ti, all the Fourier amplitudes u nm {k) and v nm {k) 
are irrelevant, they decrease when the cutoff k is moved towards zero. For f3 2 < 8tt, at least 
one of the Fourier amplitudes becomes r elevant. This phase s tructu re is very similar to 



that obtained for the massless SG model (Nandori et. al. 



2001 



2004J). Nevertheless, there 
is a significant difference. In case of the generalized LSG model, the coupling constant J k 
is always a relevant parameter and, consequently, the linearization loses its validity with 
decreasing scale k for any value of (3. Therefore, one has to determine the WH-RG flow with 
taking all the non-linear terms into account in order to obtain reliable IR behaviour for the 
coupling constants. The numerical solution of the problem is in progress. 

Finally, the UV scaling laws obtained above and applied to the LSG model (J3J) via 

uio(k) = uoi(k) = u(k), (24) 



are compared to t 



Pierson et. al. 



ie flow given by Eq. (J3J) found in Ref. ((Pierson and Valla [1,992, 



19921 ) in the dilute gas approximation. The linearized WH-RG equations 
((22J) applied to the LSG model reduce to 

p 2 



0. 



u 



An 



k^-J=-2J. 
dk 



(25) 



These turn out to be the same as those of Eq. (jlj) found in Ref. ((Pierson and Valid Il992) 
except of the loss of the scale-dependence of (3 due to the usage of the LPA. Therefore, our 
approach with the non-linear terms included in the RG equations offers the possibility to 
achieve an improvement in the determination of RG flow as compared to the results of the 
dilute gas approximation. 



VI. SUMMARY 



The new qualitative picture of the vortex length-scale dependence obtained by recent 
current transport measurements on BSCCO high temperature superconductor (HTSC) sin- 
gle crystals indicates the need for a better description of the critical behaviour of vortices in 
layered systems. In this paper we investigated the phase structure and the critical behaviour 
of the layered sine-Gordon (LSG) model which is assumed to be one of the models describing 
the vortex dominated properties of HTSC materials reasonably well. The renormalization of 
the LSG model has been considered by the differential renormalization group (RG) method 
using a sharp momentum cutoff. The usage of sharp cutoff is a necessity if a spinodal in- 
stability occurs during the renormalization of the model. The Wegner-Houghton (WH) RG 
equation as well as the tree-level blocking relation (which one should use in case of spin- 
odal instability) have been derived in the local potential approximation for the LSG model. 
The UV scaling laws and the phase structure of the LSG model have been discussed at the 



linearized level by the WH-RG method. The well-k nown UV scaling 



dilut e gas approximation previously in the literature ((Pierson and Vail 



aws obtained in the 



1992 



Pierson e,t. al 



19921 ) were recovered. Finally, it was argued that our approach offers an improvement in the 
determination of the RG flow as compared to that found in the dilute gas approximation 
when all the non-linear terms are taken into account in the WH-RG equation. 
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